About possible Phonon to Magnon alignment in 2 dimensions and theory of 
superconductivity in Copper-Oxide planes. 
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We suggest that the phonon dispersion in cuprates becomes strongly anisotropic due to interaction 
with spin waves; moreover the phonon dispersion becomes singular along \k x \ = \k y \ directions. This 
would allow more electrons to form Cooper pairs and increase temperature of the superconducting 
transition. The interaction of phonons with spin waves is more important than the interaction of 
phonons with free electrons, because spin waves do not have the Fermi surface constrain. 
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The motivation for this work is to look for phonon- 
mediated mechanism of the high-temperature supercon- 
ductivity which would give d-wave order parameter ob- 
served experimentally^^, particularly by scanning tun- 
neling imaging of non-magnetic dopants.^ This letter 
will be focused on the dispersion of longitudinal acoustic 
phonons interacting with spin waves. This interaction 
has been introduced earlier to explain Raman spectra 
of high T c cuprates^; it is similar to interaction of opti- 
cal phonons with magnons^ seen possibly in midinfrared 
absorption of cuprates 6 ' 7 , Raman spectroscopy of Spin- 
Peierls systems^, and reflectivity of CuO crystals^. We 
basically undertake a search for new anisotropic phonon 
modes seen in photoemission^fi 

The electron-phonon interaction in the Fermi-liquid 
theory does not lead to interaction between phonons and 
spin waves. The phonon propagator is renormalized by 
the polarization bubble of electron gas, the latter is pro- 
portional to the electron density in the same way as 
conductivity and diffusion constant. However electron 
preserves the spin when interacting with phonon so spin 
waves cannot interact with phonons. 

Alternative description of spin waves - magnons can 
be derived from the Heisenberg Hamiltonian, just by at- 
tributing local spins to sites. In general the interaction 
energy J ~ U 2 /t, where U is the charging energy of 
local electron and t is hopping matrix element. This 
matrix element depends exponentially on inter-atomic 
distance and so does the spin interaction J. We will 
show later that the spin-phonon interaction splits the dis- 
persion of the longitudinal two-dimensional phonons and 
new branch becomes singular and anisotropic, namely it 
goes like 



, 1 dJ y 

J da | cos(k x a) 



cos(k y a)\ 
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and the divergence is cut off by the condition u> <C 
J [sin 2 (k x a/ 2) + sin 2 (fc y a/2)], where a is the lattice con- 
stant of square lattice. The simple form Eq.^is valid for 
lj ^> u>o(k), where u>o(k) is the dispersion law of longitu- 
dinal acoustic phonons. 

The extra lattice stiffness comes out of zero-point en- 
ergy of spin waves. The effect is similar to the nature 



of attraction between bodies (Van der Waals forces) due 
to zero-point oscillations of electromagnetic wavesi^i By 
analogy, one can compute magnetostriction force due 
to polarization of magnon vacuum. However this force 
would compete with all other strains in solid state and 
would be hardly observable. In the same time dynami- 
cal lj > part of the magnon polarization operator has 
dramatic effect on phonon dispersion. The right way to 
integrate out all almost standing magnon waves aligned 
with phonon wave is to solve the Dyson equation as it is 
done below by making use of zero temperature diagrams. 

Two main ingredients are necessary to conduct the 
calculation: the spin-spin correlation function and spin- 
phonon Hamiltonian (also known as magneto-elastic 
Hamiltonian). These two are not independent, for ex- 
ample taking into account next-neighbor exchange inter- 
action would affect both spin-spin correlation and spin- 
phonon interaction. The subject is even more compli- 
cated in our "target" systems - high T c cuprates; their 
spin-spin correlation function depends on doping and 
changes across the antifcrromagnct to superconductor 
transition^. Fortunately, the spin waves polarization op- 
erator contains u> integration which completely removes 
this dependence, so our result Eq.^is robust and should 
work even if the ground state does not have magnetic 
order at zero temperature. 

According to the BCS theory! 3 , the energy gap is 
proportional to the phase volume available for scattering 
for two electrons with moments k and — k on the Fermi 
surface, when they loose energy lj. The singular phonon 
dispersion leads to dramatic consequences for electron- 
electron pairing. The half-filled tight-binding model E = 
t cos(k x a)+tcos(k y a) has rectangular Fermi surface k v = 
±Tr/a±k x . Assume that electron is on the line k x + k y = 
iv /a. Then the interaction with dispersion Eq. Q will 
keep electron on the same bond kf + kf — 7r/a, since it 
allows the process with | kf — k x | = | k\ — k v | . Twice more 
scattering out phase volume is available from hot spots 
(0,±7r/a), (±?r/a, 0). For example from the point k = 
(0, n/a) the electron can be scattered out by new phonons 
to both kf + kf = it /a and kf = kf — it /a bonds. The 
order parameter becomes anisotropic Ag gft g = 2A^^ h g . 
Phonons Eq. @J allow deviations from the exact law \ kf — 
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k x \ — \kf — k v \ and the Fermi surface is not rectangular 
for non-zero chemical potential so real life Ag is smooth 
function having maxima at hot spots. 

The alignment of phonons to the crystalline lattice is 
suggestive solution to the HTSC problem. It allows to 
interacting electrons to go up hui max out of the Fermi sur- 
face ( below h = 1), increasing both the energy gap and 
transition temperature. Here w max is maximal (cutoff) 
frequency in Eq. Q ; it sets up new energy scale between 
phonons and magnons; w max is much larger than the De- 
bye frequency lod and much lower than the spin exchange 
energy J. The accurate calculation of requires regu- 
larization of the singularity in Eq. Q and it is not yet 
done. 

Assume that the antiferromagnetic coupling constant 
J > depends on distance 

H = ^2 JSiSi' = ^2(J + - Uii)av i)SiSi> (2) 
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FIG. 1: The Dyson equation for the collective phonon - two 
magnon excitation. Vorteces have complicated k, q depen- 
dence, roughly like k 2 ' 5 . 



where each pair of adjacent sites (ii') on square lattice 
is counted once, vector a^i points from i' to i, Ui is dis- 
placement of ith site, and Si is the spin on ith site. 

Consider chain of sites counted along x axis for given 
y, then the interaction term is 



(3) 



for bond j parallel to x-axis, and y-components of k 
should be taken if the bond j is parallel to y axis, be- 
sides [b^jbt,] = Sftp are boson operators, £ = (j, t), and 



% = {sm(k x a/2),sm(k y a/2)) 



(6) 



The propagator for this real field is 



For the first glance this interaction vanishes because it 
is proportional to SW X S — V x 5 2 = 0. Extracting the 
following identity from the spin product 



[S 



i+i 



Si~i] [Si+i — Si-i] — 



we will find interaction being proportional to higher 
power of derivatives. After some algebra 
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y^(Mz - Ui')ai>i(Si> - SiY 



(4) 
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0(&O = -W(O^')> 
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U) 2 — 07q (k) + id 



(7) 

In further computations wc will not need the explicit form 
of u>o(k), and the only assumption is that < u>o(k) < 
u>d, where u>d is the Debye frequency. In other words we 
need u>o(k) to be limited for all k. 

The quantization of spin waves is convenient to start 
with choice of z axis parallel to the direction of the mag- 
netic order parameter; then the commutation relation of 
spins [S x , S y ] = i/2 leads to field operator 



that can be verified by expansion of the square term. 
The magneto-elastic energy derived here is invariant un- 
der spin rotation as opposite to the so-called relativis- 
tic termii The relativistic part of the magneto-elastic 
energy is responsible for acoustic Faradey's effect - rota- 
tion of phonon polarization around magnetization axis or 
Neel's vector, observed recently in two dimensions^ It 
will not be discussed here because i) the effect does not 
exist in disordered phase, ii) it is less important for large 
k. 

For each bond j = (ii'} we will introduce the time 
dependent phonon field 4>{j) — (Hi — i/V)aV iy/p, where 
p is the material density. For longitudinal phonons the 
field is quantized as 
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2n x k x Ik - 
v k --{b?e lkr > 



-Iujq {k)t 



c.c.} (5) 
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and the propagator of spin waves 



■iuj s (k)t 
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(8) 



lu — Lu s (k) + id 



Now we can put z axis perpendicular to the plane of 
atomic lattice, as it was before spin wave quantization. 
The magnon dispersion for the Hamiltonian JSiSy is just 
u) s (k) — 2J«;§, and we should preserve it since it should 
be consistent with interaction Hamiltonian. 

The interaction Hamiltonian is the sum over a = x,y 



k,r k u q-k/2 u ' (r+ k/2 



c.c. 



(10) 
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where 



^ fc == / = (n) 



2pw (fc) 



Then the Dyson equation for the phonon propagator be- 
comes matrix, 



1 



D % 



J a/3 u % 



(0) " «3 



and integration of the magnon loop, see Fig. ^ is our 
primary task. 

The frequency integration is simple and keeping real 
part in the answer we have 



du' 



1 



(13) 



The <f integral in turn is complicated. At low frequencies 
important for phonon dispersion 



wfe k J k ' k ^k^k 



2-nJ 2 a 2 



and it has the pole 



(14) 



sin 2 M» / S i n 2 ^ 
4 = — Tkh .-2 la + non-singular term , (15) 



where sin 2 > sin 2 and one should interchange k x 
with k y when sin 2 < sin 2 

There are two branches of the phonon dispersion, as 
obtained from Eqs. (ttlll . and detD -1 - = 0. Near the 

pole of one branch is the original phonon dispersion 
lu = ujo(k) and another branch is 



LJ 2 — 0J 2 (k) 
(%) 2 



-u,[/f +/f]/ jE + n.s.t. = 



(16) 



where n.s.t. stands for non-singular terms. Our result 
Eq. with all prefactors reads 



(-k) 2 [ir+mh 



(17) 



valid for co <C J(kj-) 2 - condition for w-expansion of M - . 

In conclusions the phonon dispersion is modified by in- 
teraction with spin waves; the effect is strong for large 
phonon wave numbers (new branch goes oc fc 3 ) and 
should be observed experimentally. The phonon disper- 
sion becomes singular along k x = ±k y lines opening a lot 
of room for electron pairing, that could be key solution for 
high-temperature superconductivity. This work has been 
inspired by discussion at Physics department of Stanford 
University. The experiment showed that low k part of 
the phonon spectrum depends strongly on electron den- 
sity where as high k part is completely insensitive. This 
is in agreement with our results, because the renormal- 
ization of D^°l by electron-phonon interaction does not 

ujk J i 

change u>(k) given by Eq. i|17fl. 
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